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Abstract 

The role of strong absorption of particles in intermidiate and final states has been 
considered. The range of applicability of phenomenological model of absorption has been 
studied. This model is nonuniversal. Its applicability depends on the type of interaction 
Hamiltonian and matrix element used. We also demonstrate that the violation of the 
unitarity condition can produce a qualitative error in the results. The absorption (decay) 
in the final state does not tend to suppress the total process probability as well as the 
probability of the channel corresponding to absorption. This is true for the reactions, 
decays and nn conversion in the medium. 
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1 Introduction 



In [1] it was shown that in field-theoretical and phenomenological models the effect of final state 
absorption acts in the opposite directions. In the following this problem is considered in detail. 
We adduce an additional arguments and study the reasons for this disagreement. This also 
makes sense if one considers that some problems were solved by means of the above-mentioned 
phenomenological models only. Also we study the range of applicability of phenomenological 
models. 

Phenomenologically, the absorption is described by an optical potential [2] . For illustration, 
let us consider a free-space decay a — > bn, for example, A — > nn°. For a decay in nuclear matter 
we have 

a^b + n^b + M (1) 

(M are the annihilation mesons) because n annihilates in a time r ~ 10~ 24 s. 

By way of another example we consider the nn transitions [3-5] in nuclear matter followed 
by annihilation 

n^n^M. (2) 

The antineutron annihilation should be described by an Hermitian Hamiltonian 7i a . In the 
phenomenological models 

H a ^H = tlmUn^n^n, (3) 

where U n is the optical potential of n, H is the phenomenological absorption (annihilation) 
Hamiltonian. For brevity, ReU n will be omitted, except if otherwise noted. 

In practice, the absorption and decay are described by a distorted wave [6], or dressed 
propagator (see, for example, refs. [7,8]). To study the model as a whole, one should write the 
total interaction Hamiltonian Tij. In specific calculations the Hamiltonian H,, as a rule, is not 
adduced. However, the corresponding terms in the distorted wave or Green function originate 
from 7i. Due to this, we consider the problem at the level of the effective (not fundamental) 
Hamiltonians. 

In the case of process (1), the phenomenological model is given by 

n I = n 1 + n a ^n 1 + n, (4) 

where Hi is the Hamiltonian of the free-space decay a — > bn. 

The phenomenological interaction Hamiltonian of process (2) is 

T~ti = Ti-nn + *H, 
Hun = e^n^n + H.c, (5) 
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e = 1/t. Here Ti, nfi is the Hamiltonian of nn conversion [5], r is the free-space nn oscillation 
time. As we will see later, process (2) is an ideal instrument for the study of the final state 
absorption and we focus on this process. 

On the one hand, model (3) is very useful because it greatly simplifies the calculation. 
On the other hand, the Hamiltonian H is non-hermitian and so model (3) is effective. Its 
applicability range is restricted. 

We consider the decay a — > bn and the nn transition in the medium and elucidate what pro- 
cesses can be described by the effective Hamiltonians (4) and (5). More complicated processes 
are considered as well. In other words, we study the range of applicability of model (3). We also 
study the suppression of the processes mentioned above due to final state absorption. This is a 
question of principal because the calculations with hermitian and non-hermitian Hamiltonians 
give opposite results. 

For our purposes it is sufficient to consider the simplest potential U n =const. We per- 
form concrete calculations and show that the unjustified application of a model can produce a 
qualitative error in the results. This is primarily true for the total probability of decays, nn 
transitions and reactions. Formally, models (4) and (5) can lead to an additional suppression 
of the total process probability as well as of the probability of the channel corresponding to 
absorption in comparison with the calculations with hermitian 7i a or, similarly, calculations 
with hermitian 7i a can tend to increase the above-mentioned values. 

With the substitution ilmU n = —iT x /2, where T x is the width of some free-space decay 
n — > x, the effective Hamiltonians (4) and (5) describe the free-space two-step processes: a — > 
b + n — > b + x and n — > n — > x, respectively. So, when referring to Hamiltonian 7i, we keep in 
mind the decay as well. 

The paper is organized as follows. In sect. 2 the simple but important statement related 
to the final state interaction is proven: the opening of a new channel leads to increase the 
total decay probability. It turns out, that for the nn conversion in the medium model (5) 
contradicts this statement (sect. 3). The same is also true for the free-space decay (sect. 4). 
Section 5 shows that the reason for this is the non-unitarity of the S'-matrix and the structure 
of the Green function. In this connection we review the origin of the complex self-energy E 
in quantum electrodynamics (QED) and optical potential theory and point out the principal 
distinctions with respect to the model under study (sects. 5 and 6). The value and physical 
meaning of ilmU for various Hi are analyzed as well. In sect. 7, we qualitatively discuss 
more complicated Hamiltonians and matrix elements. Field-theoretical and phenomenological 
approaches are compared in sect. 8. The results are summarized and discussed in sect. 9. 
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2 Absorption in the final state 



To clarify the role of final state absorption, we prove a simple mo del- independent statement. 
We consider the decay a — > bn in the medium. Let H s and H a be the hermitian Hamiltonians 
of the scattering and annihilation of n, respectively. The total a-particle decay probability 
W t is W t = W n + W a , where W n and W a are the probabilities of finding an antineutron and 
annihilation mesons M, respectively (see fig. 1). 




Figure 1: (a) The decay a — > b + n in the medium without considering annihilation, (b) The 
same as (a), but with annihilation. The antineutron annihilation is illustrated by a circle. 

Let H a = and H s ^ 0, then W t = W n (see fig. la). Now, let us turn on the perturbation 
H a (see fig. lb). In the lowest order in H a we have 

W t (H 8 + H a ) = W n (H s + H a ) + W a (H s + H a ) > W n (H s + H a ) = W n (H s + 0) = W t (H s + 0). 
Thus, 

W t (H s + H a )>W t (H s + 0). (6) 

In the equality W n (H s + H a ) = Wn(H s + 0) it was taken into account that only terms of zero 
order in H a give a contribution to W n . Inequality (6) can be written in terms of decay widths. 
We use the probabilities for reasons given in sect. 5. 

Similar process for the nn transition [3-5] in the medium is shown in fig. 2. Obviously, for 
this process, inequality (6) is true as well. 

In eq. (6) we can put H s = 0. Instead of the annihilation we can consider any process, 
for example a decay in the final state with the Hamiltonian H x . It is important that there is 
no interference between diagrams a and b, and H x is small. The simplest case is given by fig. 
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Figure 2: The same as fig. 1 for the nn conversion in the medium. 

1, where 7i s = and the circle corresponds to the decay. Then in the lowest order in 7i x we 
have W t (Hi + H x ) > W t (Hi + 0). The expression H s + H a is nothing more than a notation. 
It illustrates the presence of two channels. In principle, the absorption and scattering can be 
described by one and the same interaction Hamiltonian as in QED. 

Inequality (6) shows an obvious fact: the opening of a new channel (annihilation) leads to 
increase W t . Obviously, this is generalized to more complicated processes: reactions, decays 
and ab conversion involving final state absorption. 

It turns out that models (4) and (5) give the opposite result. This can be easily shown for 
nn transitions in the medium. 

3 Absorption in the phenomenological model 

In the standard approach (see, for example, refs. [9-13]) the nn transitions in the medium are 
described by Schrodinger equations: 

(id t — H )n(x) = en(x), 
(idt — H — V)n(x) = en(x), 
H = -V 2 /2m + U n , 
V = U n -U n = ReU n - iT/2 - U n , (7) 

n(0, x) = 0. Here U n and U n are the potentials of n and n, respectively; e is a small parameter, 
T being the annihilation width of n. 

For V =const. in the lowest order in e the overall nn transition probability (the probability 
of finding an n or annihilation products) in a time t is [13] 

W t {t) = 1- | U u {t) \ 2 = 2ImT«(t)- \ T u {t) 2ImT ii (i), (8) 
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T«(t) = i{t/Vf[l - iVt - exp(-iVt)], 

where U(t) is the evolution operator; Uu(i) = 1 + iT u (t) =<n(0) \n(t) >. 

If V =const., system (7) has an exact solution. Since e is extremely small, only lowest order 
in e is commonly taken into account. This is a sole approximation made in the calculation of 
the Tu in the framework of model (5). 

At least for small V 

W t (ReV + ilmV) < W t (ReV + 0), 

dwydr < o, (9) 

which contradicts to (6). Indeed, let Tt^> I. Then 

ty ' (i) - 2e VvT + (r/2)^ 4A/r - (10) 

This is a well-known result [5,11,12] for nn transitions in nuclear matter. If (T/2) 2 > (ReV) 2 
(the realistic set of parameters fits this requirement), dW t /dT < 0. At the point ReV = 0, 
dWydr < as well. 

In the opposite limiting case | Vi|-4C 1, 

W t (t) =e 2 t 2 (l-rt/6) (11) 

and we arrive at eqs. (9) again. On the other hand, at small V inequality (6) is also valid. 
Thus (11) contradicts to (6). In model (5) the effect of absorption acts in the opposite (wrong) 
direction, which tends to the additional suppression of the nn transition. 

In (6) and (9) physically identical procedures have been done: 7i a = — > 7i a ^ and 
\mV = — > ImV 7^ 0, respectively. The results are opposite. Equation (10) shows that the 
potential ReV suppresses the nn transition, which is certainly correct, however, T acts in the 
same direction, which seems wrong. 

To clarify the structure of (10), we consider the same problem by means of a diagram 
technique [1] (see fig. 3a). Here we use the S- matrix rather then an evolution operator. Put 
U n = ReU n = for simplicity. The Hamiltonian (5) has the form 

Hj^fA + if.c-^fA. (12) 
The antineutron propagator G and total process probability Wf(t) are 

G = l/(e n - p 2 /2m + iT/2) = 2/iT, 

W t d (t) = -2lmeGet = (13) 
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Figure 3: (a) The on-diagonal matrix element Ta corresponding to the nn transition in the 
medium. The Hamiltonian of the n-medium interaction is given by eq. (3). (b) The same as 
(a) for the decay a — > bn. 

where p = (e n ,p) is the neutron 4-momentum; e n = p 2 /2m. So Wf = Wt, where Wt is given 
by (10). The T-dependence of W t is conditioned by the propagator G. This fact is common for 
a 2-tail diagram (fig. 3a) and diagrams with a momentum transferred q ^ (figs, lb, 3b and 
4b). 

In this section the standard scheme of calculation has been used. It is based on the Hamil- 
tonian (5) and equation W t = 2\mTu. This is the sole way of calculation of W t in a one-particle 
model. For brevity, this model will be denoted as model (5). We, thus, see that (5) realized by 
means of the equations of motion or diagram technique contradicts to inequality (6). 



4 Free-space process 

To avoid questions connected with the medium corrections, we consider the imaginary free-space 
process 

n — > n — > pe + u, (14) 

in which the neutron decay is excluded. The hermitian Hamiltonian is Hi = e^n^n + H.c.+H^, 
where 7if is the Hamiltonian of the free-space /9 + -decay n — > pe + u. The corresponding diagram 
is shown in fig. 4a. 

Alternatively, if we want to use model (5), we have 

Hi = + H.c. - i^A, (15) 

where is the width of the free-space /9 + -decay n — > pe + u. Comparing with (12), it is seen 
that we can use all the formulas given above in which r = Tp. The total probability of the 
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a b 

Figure 4: (a) The free-space process n — > n — > pe + z/. (b) The diagram corresponding to the 
effective amplitude of the process (14) (see the text). 

free-space nn transition Wf is given by (10) or (13): 

Wf(t) = -2lmeGet « 4e 2 */I>, (16) 

r^t ^> 1. The free-space nn conversion is drastically suppressed by the decay in the final 
state. Indeed, the free-space nn transition probability Wf is Wf = e 2 t 2 (see (11), where T = 0) 
and correspondingly Wf(t)/Wf(t) ~ 1/r^ < 1. This is clearly wrong because the state the of 
intermediate n (see fig. 4a) coincides with the final state of the free-space n — > n transition, and 
so the /3 + -decay makes no influence on the subprocess of the nn conversion. This is sufficient 
to reject the model (5). 

We should make a small comment. The process shown in fig. 4a represents two consecutive 
free-space subprocesses. The speed and probability of the overall process are defined by those 
of the slower subprocess. Since l/Tp <C t, the /5 + -decay can be considered instantaneous: for 
any t± < t the (3 + -decay probability Wp is Wp{t\,t) « I. Then, the total process probability 
W 13 is defined by the speed of the nn conversion: W p « W f ~ t 2 instead of ~ t/I> 

Let us try to compose an effective model which produces Wf through the direct calculation 
of the off-diagonal matrix element. We consider fig. 4b. It differs from fig. 4a by the fact that 
the antineutron is in the potential V = —iTp/2. (Nonsense, of course.) The process amplitude 
is M e ff = —eGMp, G = — 1/V, where Mp is the amplitude of the /5 + -decay. For the process 
width we have T e ff — J d$ \ M e ff \ 2 /2m = 4e 2 /r /3 , which coincides with Wf/t. So we have 
fixed <f \ =<pe + v\ and found that the effective amplitude of process (14) which produces the 
result (16), is given by fig. 4b. 

The fallacy of this model is obvious. Certainly, this is an illustration only, but structure 
(16) can be obtained by means of the Green function G = —1/V ~ 1/T solely. 

In the calculation of the W t ,Wf and Wf model (5) is used. The result W t ~ 1/r is very 
sensitive to T. However, the T-dependence of Wt contradicts to inequality (6). Besides, result 
(16) is unrealistic. Therefore, this model should be revised. 
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5 Unitarity and self-energy 



In this and next sections the reasons for disagreement indicated above are studied. If T = 0, 
system (7) is certainly correct. Consequently, it is necessary to revise the role of ilmUn- This 
question has been considered in [1]. Taking into account the importance of this problem, we 
adduce more direct evidence using the {/(t)-operator only. The approach based on the evolution 
operator is more general than the ^-matrix one, since in this case the time- dependence of the 
process does not need to be Wt = 1 — exp(— T t t) (see (8)). Also it is infrared-free, which is 
essential for ab transitions [13]. 

The non-hermit icity of Ti implies that 

(U(t)U + (t)) f i = 6fi + a fi (t), (17) 

otfi 7^ 0, resulting in 

W t (t) = J2 I T fi(t) I 2 - 2ImT ii (i) + a u {t) ^ 2lmT u (t) 
because the value of 2ImTjj is extremely small: 

2ImT^ ) = ^ < 10- 31 , 

where the standard set [1, 14-16] of parameters e, t and T has been used 
(8) is invalid. 

For the S-matrix, the conclusions are the same: (a) The basic relation 

]T | T fi 2ImT ii (20) 

is inapplicable, (b) The physical meaning of ImS = — T/2 is uncertain because it is clarified 
using relation (20). We would like to emphasize this fact. 

On the one hand, the nn transition probability is very small (see (19)), and on the other 
hand, the term ilmUn plays a crucial role because it enters the leading diagram (see (13)). 
Because of this for the problem under study the unitarity of the ^-matrix is of particular 
importance. 

Thus, the non-hermitian Hamiltonian (3) leads to inverse T-dependence of W t and to the 
imaginary self-energy. In QED the Green function above threshold contains an imaginary self- 
energy ImS 7^ as well. However, in the case of QED the situation differs principally. ImS is a 
complicated function of parameters of the hermitian Hamiltonian. It appears at higher orders 
in a. The width T makes its appearance after a Dyson summation of the relevant self-energies. 
In order to correctly enforce unitarity, the notation of the "running width" has been introduced. 



(18) 



(19) 

. We thus see that 
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The importance of unitarity condition is well known [17,18]. Nevertheless, the non-hermitian 
models (3)-(5) are frequently used for the reasons given in sect. 1. In particular, all existing 
calculations of nn transitions in the medium are based on model (5) (see, for example, [13] for 
future references). 

With the substitution ilmU n = —iT x /2, where T x is a width of some free-space decay 
n — > x, the Green function (13) describes the non-relativistic resonance; Hamiltonians (4) and 
(5) correspond to the free-space two-step processes: a — >b + n^b + x and n — > n — > x, 
respectively. This is obvious because the absorption can be considered as the decay of a one- 
particle state. Formally, in these cases all the results are also true. Nevertheless, the resonances 
invite an additional consideration. As far back as 1959, M. Levy remarked that there does not 
exist a rigorous theory to which various phenomenological methods of treating resonances and 
decays can be considered as approximations [19]. Attempts have been made at an axiomatic 
theory [20,21]. 

The above-mentioned difficulties take place for absorption as well. These conceptual prob- 
lems are beyond the scope of this paper. We deal with concrete models (3)- (5) and hence 
propagator (13) because they are frequently used. As for resonances and decays, we only draw 
the formal analogy between absorption Hamiltonian (3) and phenomenological Hamiltonian of 
decay -iY/2^. 

We also note that decay (14) can be calculated by means of the usual field-theoretical 
approach, but the problem should be formulated on the finite time interval [22] since fig. 4a 
contains an infrared singularity. 



6 Optical potential 

The problem is not only in the unitarity. It is in the correct description of the absorption on 
the whole. In the theory of optical potential ilmU is non-hermitian as well. However, the 
picture differs principally in this case. In this section we compare the equation of motion and 
the problem under study from the standpoint of the use of an optical potential. 
In the case of Schrodinger equation 

(id t -H - U n )n = (21) 

the scheme is as follows. Since U n is non-hermitian, the condition of probability conservation 

1 =| U« | 2 +W Sch (22) 

is imposed. Here Ws c h is the loss of n intensity. The matrix element of evolution operator is 
found to be 

U u (t) =< n(0) | fi(t) >= e - iReC/ «* e ImC7s *. (23) 
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From (22) and (23) it is seen that: (a) W Sch 7^ if and only if ImU n < (or T > 0, when 
lmU n = — r/2). (b) If T increases, W^di increases as well: 

^ = ^(1 _,-»)> a (24) 

This agrees with (6) and is in contradiction with (10), (11) and (16). 

The procedure given above is based on two points: 1) In (21) \mU n has a clear physical 
meaning. It is defined by the continuity equation corresponding to (21). 2) The additional 
bound (22) provides the probability conservation (unitarization). By means of (21) and (22) 
U n is fitted to the p-atom and low energy scattering data. 

For more complex problems these requirements, as a rule, are not fulfilled. We demonstrate 
this for model (5). The fit of (7) and (8) is impossible since there are no experimental data. As 
a result we have (18) with the consequences considered above. In addition, we try to realize 
the scheme given for (21). 

The coupled eqs. (7) give rise to the following equation 

{8 2 t + id t (V + 2H ) - H 2 - H V + e 2 )n(x) = 0. (25) 

According to (8), n(x) is sufficient to get W t . 

Even the first step of the scheme described above is not realized: one cannot get the conti- 
nuity equation from (25). The S- matrix consideration accomplishes nothing because eq. (20) 
is inapplicable. 

Equations (7), i.e. model (5), describe only W n . In this case U n can be included in the 
distorted wave of the antineutron which is the eigenfunction of eq. (21), and this justifies the 
model. 



7 Generalization 

If instead of Hamiltonian (5) we take 

Hi = U r , d + H, (26) 

where 7i r and Ttd correspond to free-space reaction and decay, respectively, the qualitative 
conclusions do not change because the heart of the problem is in the Hamiltonian H. As an 
example, let us consider the decay a — > bn in the medium. Let T n and T a be the widths of 
decays with n and the annihilation mesons in the final state, respectively; r t is the total decay 
width, F t = T n + r a . The corresponding partial decay probabilities are Wn A ~ T n a t. W n and 
W a are the same as in sect. 2. To draw the analogy to nn transitions, we use the probabilities 
W. 
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Equations (7) are time-dependent and so the evolution operator has been applied. For the 
decays the ^-matrix is used. In (18) one should replace T(t) — > T. The interaction Hamiltonian 
is given by (4). We have 

r t = (27) 

where To is the normalization time, To — > oo. The matrix element T„ is shown in fig. 3b. In 
principle, the antineutron propagator in the loop should be calculated through the hermitian 
Hamiltonian H a : G = G(H a ). Model (4) means that 

G(H a ) - G(H) = G(-<r/2) = (28) 

Pn - m + li /2 

where p n is the antineutron 4-momentum. Obviously, for the matrix element shown in fig. 3b 
eq. (18) takes place as well. Relation (20) is invalid; the physical meaning of ImS = — T/2 is 
uncertain. 

The probability of finding an antineutron W n is described by an off-diagonal matrix element. 
In the distorted wave impulse approximation the interaction responsible for the absorption is 
included in the antineutron wave function: 

fi(x) = fi-l/2 e -i[(p 2 /2 m +ilml/)^px] > ( 2Q ) 

The corresponding diagram is shown in fig. la, where the antineutron state is described by 
(29). The wave function n(x) is the eigenfunction of eq. (21), which justifies the using of model 
(3) in the calculation of n(x) and W n - 

The probability of finding the annihilation products is obtained from 

W a = W t - W- n . (30) 

Since eqs. (20) and (27) are inapplicable, W t and W a are uncertain. 

We thus see that (4) describes only W n . The result is the same as for the nn transitions 
considered above. Obviously, in the strong absorption region W n -C W a and W n <C W t . 

Figures lb and 2b correspond to absorption in the final state. Model (3) is also used for the 
description of the absorption in the intermediate state (see fig. 5). The interaction Hamiltonian 
of the process shown in fig. 5 has the form 

Hi = H x + H + U{. (31) 

The quantitative study of models (26) and (31) is subject of a separate investigation. Here 
we consider only a qualitative picture. The amplitude corresponding to fig. 5 is given by 

T 5 = -M x G{iV)Mp, 
G(iT) = [(p - g - m) - (p - q) 2 /2m + ^r/2]- 1 . (32) 
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Figure 5: The decay a^b + n^b + p + e + + v in the medium. The bold line signifies the 
antineutron annihilation in the intermediate state. 

Here M 1 and Mp are the amplitudes of decay (1) and ,3 + -decay, respectively; p and q are the 
4-momenta of particles a and b, respectively; m is the antineutron mass. 

The antineutron propagator G(Ti. a ), calculated through the hermitian Hamiltonian Tt a , 
contains the loops. This leads to suppress the amplitude T 5 . G(iT) from (32) acts in the same 
direction: the probability of finding the /9 + -decay products is W 5 ~ Tp/T 2 . W 5 is less sensitive 
to model (3) than W t because the unitarity condition is not used. In this case at least there 
is no qualitative contradiction. (This question needs quantitative consideration. Most likely 
model (32) yields too great suppression.) 

Below we consider the most realistic case T 3> Tp. In the lowest order in the probability 
of finding an antineutron W n is the same as for Hamiltonian (4). For the on-diagonal matrix 
element T u and total decay probability W t the calculation scheme and conclusions are also 
identical to those for (4) and (26). The fact that the antineutron propagator in the loop is 
defined by 7if and H is not principal because the heart of the problem is in the ilmUn- 

Similarly to (30) we have 



Since eq. (20) is inapplicable, W t and W a are uncertain. 

For model (31) we conclude: (a) W n should be described correctly, (b) The major decay 
characteristics W t and W a are not described, (c) For the process shown in fig. 5 model (32) 
can be used as a first approximation. 



In sects. 2-4 on the basis of a general reasoning we concluded that in the phenomenological 
model the T-dependence of W t is wrong. Below we consider the unitary model and calculate 
directly the off-diagonal matrix element by means of the diagram technique. 



w a = w t -w n - w 5 . 



(33) 



8 Unitary model 
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If rt ^> 1, the probability of finding an annihilation mesons W a is much greater than 
W n . However, the phenomenological model describes W, n only. Recall that for the total nn 
transition probability the phenomenological model gives W t ~ 1/r (see (13)). Since W n <C W t , 
W a depends inversely on Y as well: 

W a = W t - W n « W t ~ 1/r. (34) 

For the processes which are described by Hamiltonians (26) and (31) it is sufficient to recall 
that the W t and W a are uncertain for the reasons given above. In our opinion, with correct 
consideration of the corresponding loops we will obtain dW t /dT < 0, as with (34). 

The direct calculation of off-diagonal matrix element gives the inverse T-dependence dW/ dT > 
0. Indeed, we consider the process (1). The a-particle and n are assumed non-relativistic. The 
wave function of the o-particle is = (2go^) -1 ^ 2 exp(— iqx), where q is the 4-momentum 

of the particle. As with Ti, nn , the decay Hamiltonian is taken in the scalar form Ti, 1 = 
e'^ n ^l^ a + H.c. and correspondingly 

Hi = e^ n ^ a + H.c.+H a , (35) 

e' is dimensionless. 

The process amplitude is given by 

M 1 = -e'- — - — - M a . (36) 

{Po -q -m)-(p- q) 2 /2m + tO 

Here M a is the annihilation amplitude, m is the antineutron mass, p is the 4-momentum of the 
a-particle. 

For simplicity assume that m^jm <^ 1, where mj, is the mass of the 6-particle. It is easy to 
estimate the width of decay (1): 

r x w e' 2 r/(27r 2 ). (37) 
The corresponding decay probability is proportional to T: 

W^ = Tit ~ T. (38) 

The index h signifies that the hermitian Hamiltonian is used. 
The width of process (2) is also linear in T [1]: 

r 2 ~ r. (39) 

For the Hamiltonians containing three terms the T-dependence of W^ is the same. Thus, 

w h a ~ w t h ~ r, (40) 
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where Wj 1 is the total process probability. 

From eqs. (40) and (34), we see that the unitary and non-unitary models lead to inverse 
T-dependence of the results. Because of this the calculations with the hermitian 7i a can tend 
to increase W a and W t . (See also eq. (12) of ref. [1].) 

9 Summary and conclusion 

We list the consequences of an unjustified use of the model based on eqs. (20) and (3). 

1) Equation (9) contradicts to (6) and (24). 

2) Result (16) is unrealistic. 

3) W t ~ 1/r, whereas W t h ~ T (see (34) and (40)). 

4) The physical meaning and value of the ImS are uncertain (see text below (20)). 
Model (3) was adapted to quite definite problems. It is justified for the problems described 

by Schrodinger type equations. It also describes the complicated processes (reactions, decays 
and nn transitions) with n in the final state. (More formally, model (3) can be applied to the 
calculation of W n corresponding to the Hamiltonians containing several terms (eqs. (26) and 

(31) , for example).) As a first approximation, it can be used in the calculation of the diagrams 
like that shown in fig. 5 with n in the intermediate state. In these cases W are calculated 
directly without the use of the unitarity condition and the calculation of T^. 

In other cases, when the interaction Hamiltonian contains several terms and the unitarity 
condition is used (eqs. (8), (13), (27) and (33), for example), model (3) is inapplicable. The 
calculation of the total process probability W t (and thus W a ) corresponding to inclusive reaction, 
decay or ab transition is impossible. The physical meaning of ImS is uncertain. The effect 
of absorption, as a rule, acts in the opposite (wrong) direction, which leads to additional 
suppression. In particular, model (5) gives rise to the dramatic suppression of nn transitions 
due to the annihilation in the final state, which is wrong. 

This paper also demonstrates the importance of the unitarity condition for any model of E 
[23,24]. The model should by unitary or unitarized. 

Finally, we touch upon the result sensitivity to model (3). It is seen from the Green function 

(32) . The T-dependence is masked by q. If q — > (2-tail) and m a = m, the problem is extremely 
sensitive to T: T 5 ~ 1/r. Alternatively, in the phenomenological model the nn conversion is 
described by system (7) which has an exact solution. For these reasons the nn transitions in 
the medium are the ideal instrument for the study of the final state absorption. 

We also emphasize the following: the absorption (decay) in the final state (figs, lb, 2b and 
4a, for example) does not lead to suppress the total process probability as well as the probability 
of the channel corresponding to absorption, in contrast to the phenomenological model results. 
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Therefore, the calculations based on unitary models can tend to increase the above-mentioned 
values. 

The author is grateful to Prof. E. Oset for helpful comments. 
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